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SUMMARY 


The available information on the spectrum of atmospheric turbulence 
is briefly reviewed. On the basis of these results, a method is developed 
for the conversion of available gust statistics normally given in terms 
of counts of gust peaks into a form appropriate for use in spectral calcu- 
lations . The fundamental quantity for this purpose appears to be the prob- 
ability distribution of the root -mean-square gust velocity. Estimates 
of the variation of this distribution with altitude and weather condition 
are also derived from available gust statistics . 


INTRODUCTION 


During the last 20 years, a considerable body of statistical data 
on the gust and load experience of airplanes in transport operations 
and in special flight tests has been collected by the National Advisory 
Committee for Aeronautics. These data have in the past, for simplicity, 
been reduced on the basis of simple airplane gust-response theory and 
the so-called "effective gust velocities" determined. In the last few 
years, it has been considered desirable to redefine these gust velocities 
in order to account more fully for the variations with altitude in the 
airplane response to gusts. The term "derived gust velocities" (ref. 1) 
has been applied to these redefined values, and the conversion of the gust 
data into this form is now largely completed (see, for example, refs. 2, 

3, and A) . The gust data in this form have been analyzed in detail, and 
the airplane gust experience for various types of operations and for 
various flight altitudes and weather conditions has been established 
(ref. 3). These results are finding useful application in the calcula- 
tion of gust loads for arbitrary flight plans. 

Even more recently, it has been desirable in many gust-response 
problems both to account explicitly for the continuous character of 
atmospheric turbulence and to account in greater detail for the airplane 
dynamics. For this purpose, the techniques of generalized harmonic 
analysis or spectrum analysis appear particularly appropriate and they 


2 


NACA TN 55^0 


are finding increasing application (refs. 6 to 15) • The application of 
these techniques to the gust-load design problem requires detailed infor- 
mation on the spectral characteristics of atmospheric turbulence. The 
purpose of this paper is to describe the information available on the 
spectral characteristics of atmospheric turbulence and, in particular, to 
adapt the large amount of statistical data on gust loads already available 
into a form suitable for use in spectral calculations. First, the avail- 
able measurements of the spectrum of atmospheric turbulence are reviewed. 
Then, a technique for the reevaluation of the gust statistics is outlined 
and some results obtained by the application of this technique are 
presented. 


SPECTRAL CHARACTERISTICS OF ATMOSPHERIC TURBULENCE 


Figure 1 contains a summary of most of the available airplane meas- 
urements of the power spectrum of atmospheric turbulence . The first of 
these measurements was made by Clementson at the Massachusetts Institute 
of Technology (ref. 6) and subsequent measurements were made by the National 
Advisory Committee for Aeronautics (refs. 10 and 11), the Douglas Aircraft 
Company, Inc. (ref. 12), and Summers at the Massachusetts Institute of 
Technology (ref. 15). The curves shown represent the various power spectra 
and were obtained under different weather conditions. The abscissa is the 
frequency argument 0 . which has the dimensions of radians per foot and is 
equal to 2it divided by A, the gust wavelength. (The data shown cover a 
range of gust wavelengths from about 10 feet to 5^000 feet.) The ordinate 
is the power density normalized to the power-density values at 0 , = 0.01. 
This normalized ordinate is used in order to facilitate comparison of the 
various spectral shapes . 

The spectra in all but one case are for the vertical component of the 
turbulence. In one case marked by the letter H, the spectrum is for the 
horizontal or longitudinal component of turbulence . Examination of these 
results indicates that the spectral shapes appear to be relatively con- 
sistent; in all cases, the power decreases rapidly with increasing fre- 
quency. In fact, in most cases, the spectra appear to be inversely pro- 
portional to the square of the frequency. This spectral shape of 1 /^^ 
is in reasonable agreement with theoretical results obtained for the 
spectral shape at the higher frequencies in the theory of isotropic tur- 
bulence. At the lower frequencies, the situation is not so clear, few 
measurements being available for frequencies of 0 , < O.OO5. Some addi- 
tional measurements obtained at the Cornell Aeronautical Laboratory do 
cover these lower frequencies and indicate a flattening of the spectrum 
in this frequency region. 

In addition to these variations in spectral shape, the various meas- 
urements also differ in turbulence intensity. The individual root-mean- 
square values are estimated to vary from roughly I.5 to perhaps 8 feet per 
second, which, as will be seen, represent the relatively lighter levels of 
atmospheric turbulence . 
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Another soiirce of information on the spectral characteristics of 
atmospheric turbulence is the measurements at lower altitudes made from 
meteorological towers. A large number of such spectral measurements 
have now been obtained. A few representative measurements from refer- 
ence l6 are shown in figure 2 and were obtained at an elevation of about 
300 feet and for various conditions of average wind speed V as shown 
in figure 2. These measurements extend to lower frequencies (longer gust 
wave lengths ) than the available airplane measurements . At the higher 
frequencies these results approximate the same l/^2^ fomm, characteristic 
of the airplane measurements. In addition, at the lower frequencies, a 
definite tendency toward a flattening of the spectra can be noted. The 
variation in the spectrum intensity with wind speed should also be noted. 

Because of the general characteristics of these spectral measure- 
ments, it has been convenient in theoretical studies to use the following 
analytical expression for the turbulence spectrum 




1 + 3U^L^ 

(l + 


( 1 ) 


where 0, - 2jt/A and A is the gust wavelength. This expression has been 

useful in wind-tunnel studies of turbulence and, as will be seen, has the 

general characteristics of the measured spectra of atmospheric turbulence . 

The equation has two parameters: the mean-square gust velocity which 

describes the overall intensity, and the so-called scale of turbulence L 
which, in a rough sense, can be considered to be proportional to the aver- 
age eddy size . Curves for this expression are shown in figure 3 foi* 
values of L of 200, 600, and 1,000 feet. At higher frequencies, these 

curves all approach a shape of l/n^ but differ in the frequency at which 

the flattening occurs. For increasing values of L, the curves flatten 
out at lower frequencies. Comparison of these curves with those of fig- 
ures 1 and 2 and other measurements of the spectrum for atmospheric tur- 
bulence has suggested that representative values of L for atmospheric 
turbulence are from at least several hundred feet to over 1,000 feet. 


For design purposes, the overall gust experience in operations is 
of concern. Presumably, the overall experience consists of various 
exposure times to each of the spectra shown in figures 1 and 2 as well 
as to other spectra associated with different weather conditions. It 
thus appears important to determine the proportion of flight time spent 
under these various conditions of turbulence. In the remainder of this 
paper, some results derived from available gust statistics on this prob- 
lem will be described. The basic approach used is described in detail 
in reference lA and will only be outlined for present purposes. 
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RELATIONS BETWEEN PEAK COUNTS AND SPECTRA 


It will be recalled that available gust-load statistical data are 
normally given in the form of counts of acceleration peaks exceeding 
given values. It is therefore necessary to relate such peak counts to 
the spectra of the turbulence. For this purpose, it appears necessary 
and reasonable to make some simplifying assumptions. First, it will be 
assumed that the turbulence for a given weather condition is a Gaussian 
process. This assumption implies that the turbulent velocity fluctua- 
tions and also the airplane response for linear systems have a normal 
or Gaussian probability distribution; this distribution appears from 
available measurements to be a reasonable approximation. The Gaussian 
character for the turbulent velocity fluctuations for a given weather 
condition will be used as a building block to construct the operational 
load history which covers many weather conditions and, in the overall, 
is far from a simple Gaussian process. Secondly, it appears necessary 
to make some simplifying assumptions on the form of the power spectrum 
of the gust velocity. On the basis of the available measurements, it 
appears reasonable to assume that the spectral form is given by equa- 
tion (l)« These assumptions and available relations between peak counts 
and spectra for the Gaussian case provide the basis for the analysis to 
be presented. 


Homogeneous Case 

For a Gaussian process y(t), the number of peak (maxima) values N(y) 
per second exceeding given values (for large values of y) is given by 
the following relation (ref . 17 ) : 



-y2/2cy2 


( 2 ) 


where 

y a random variable 

on frequency, radlans/sec 

0(cd) power spectrum of y 

a root-mean-square value of y 

When this relation is applied to the acceleration history of an airplane 
in homogeneous turbulence, the spectrum of normal acceleration for a 
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linear system is^ in turn^ related to the spectrimi of atmospheric tur 
bulence by the relation 




O) 


where 1a(“) and Oy(oi) are the power spectra of the normal accelera- 
tion and vertical gust velocity, respectively, and T(cii) is the ampli- 
tude of the airplane acceleration response to unit sinusoidal gusts 
encountered at frequency o). Equations (2) and (3) thus relate the 
mmiber of peak accelerations per second exceeding given values to the 
spectrum of the turbulence. When these relations are applied to an air. 
plane in flight through turbulence of a given spectrum (which in the 
present application can be considered to be given by equation (l)), the 
number of peak accelerations per second exceeding given values of accel. 
eration a, as illustrated by the sketch on the upper left of figure 4, 
may be expressed as 


N(a) (^) 

where 

-1/2 

cjq) T^( Oi) dm 

; (4a) 

$y(m)T2(oj)da> 


If the expression for the turbulence spectrum given by equation (l) is 
applied in equation (4a), the quantity Nq may be seen to depend upon 

the value of the scale of turbulence L and the airplane transfer func- 
tion but is readily seen to be independent of the turbulence intensity or 
root -me an -square gust velocity. Thus, for a given airplane and a given 
value of L, is a constant, independent of turbulence intensity. The 

quantity Nq has the dimensions of a frequency and can be considered to 

be a characteristic frequency of the airplane response to the given tur- 
bulence . In addition, the number of peaks N(a) is seen from equa- 
tion (4) to depend also on the acceleration level a and the root-mean- 
square acceleration Oq^ which, as will be seen, depends directly on the 

turbulence intensity. 
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Composite Case 


So far, the turbulence has been assumed to be homogeneous and 
Gaussian. Now, if the ai3rplane is assumed to fly for given flight 
times ti, t2, and tj through turbulence of different intensities but 
of the same spectral form, the acceleration histories for each fli^t 
time would differ in intensity or in root-mean-sq.uare acceleration value 
as indicated by the sketches on the lower left of figure Equation (4) 
may be extended for this composite case in order to obtain the overall 
average number of peaks exceeding given values as follows: 


-a^/ 2 cT; 


N(a) = No 


ai 


-a2/2a; 


Q'2 


% 


+ — e 
t^ 


■^3 

+ e 


-a2/2cra^^ 


( 5 ) 


where 


tijj “ t2^ + tg ‘J* t^ 


Equation ( 5 ) is essentially a sum of terms of 

"^1 

weighted by the relative exposure times — , 

% 


the form 
^2 

— , and 

tip 


of equation 
t^ 


(^) 


Continuous Case 

If the airplane is now assumed to encounter turbulence of all inten- 
sities (continuous variations in the root-mean-square gust velocity) but 
of the same spectral shape (fixed value of L for present applications), 
equation (5) may be extended for the continuous case to yield the fol- 
lowing relation: 

N(a) = Nq r f(da)e"^ ^.c^a (^) 

^0 


In this case the sum of the terms of equation (5) becomes a continuous 
integral. The expression for the continuo-us case contains the same char- 
acteristic frequency Nq, which is independent of the turbulence intensity 
the same exponential term as the earlier expressions but, in addition, 
contains the function f(cfa). This function is the counterpart of the 


■^1 





ratios 


and 


in equation ( 5 ) is the probability 
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distribution of da which can be considered to define the proportion of 
total flight time spent at values of da between ^a and da + dda. 

A sketch of f(da) is shown in figure 5- The total area under the curve 
is, of course, equal to one with the percentage of the total time between 
any two values given by the area under the curve between the two values 
as indicated by the hatched area. 

Equation (6) relates the probability distribution of the root-mean- 
square acceleration to the overall number of peak accelerations per sec- 
ond. Since gust-load statistical data are usually summarized in the form 
of number of peak accelerations per unit time (or distance) exceeding 
given values, this relation may be used to determine the appropriate dis- 
tribution of the root-mean-square acceleration f(da) from the available 
data without recourse to the reevaluation of the original records. 

For a given airplane, the root -me an- square acceleration is also 
simply related to the root-mean-square gust velocity. This relation may 
be expressed simply as 


da = AOjj (T) 


where A may be considered a transfer function between root -mean-square 
gust velocity and root-mean-square acceleration and for a given case is 
a constant depending upon the gust spectrum and the airplane response _ 
characteristics. For the gust spectrum used in the present analysis, A 
is given by the following Integral 


1/2 

( 8 ) 

The value of A for this case is seen to depend only upon the value 
of L and the airplane transfer function. Thus, the root -me an- square 
values of acceleration for a given airplane and a given value of L 
may be used directly to determine the root-mean-square gust velocities. 

Equation ( 7 ) may be applied directly to determine the distribution 
of root -mean-square gust velocity for a given distribution of root -mean- 
square acceleration. The relation between these distributions is given 
■by 


A = 


0(1 + 


T^(n)dfl 
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f(0u) = Af(cfa) (9) 

where ha = Ahy. Thus, the final result desired, the distribution of 

the root-mean-square gust velocity, may be obtained by the solution of 
equation (6) for f(ha) and the application of equation (9). 

These relations can, of course, also be used to calculate the number 
of peak loads for a given operation, ^if the probability distribution of 
the root -me an- square gust velocity fCdy) is known. In terms of the 

distribution of root -mean-square gust velocity, the number of peak 
accelerations per second exceeding given values is given by 


N(a) = No f da„ (lO) 

where Nq is given by equation (4a). Thus, in this approach, the prob 
ability distribution of the root -me an- square gust velocity is the fun- 
damental quantity and takes the place in spectral calculations of the 
counts of gust peaks used in the discrete gust calculations. 


SOME RESULTS OBTAINED IN THE REEVALUATION OF GUST DATA 


In this section, the results obtained by the application of the 
relations derived in the preceding section to available statistical data 
on gust loads are presented. Three principal problems are encountered 
in such applications and involve the choice of an appropriate value of 
the scale of turbulence L, the determination of the airplane accelera- 
tion_transfer function for gusts T(cd) (or, more particularly, the value 
of a), and the determination of the value of the characteristic fre- 
quency Nq. The value of L chosen for present purposes was 1,000 feet 

which, on the basis of available data, appears to be a reasonable estimate 
for an average value for atmospheric turbulence. In the determination > 
of the airplane transfer function and the values of A for the various 
airplanes involved, the airplane accelerations were assumed to be due 
entirely to the vertical gusts and the airplane was assumed to respond 
in vertical motion only (no pitch, rigid-body condition). A justifica- 
tion for this restriction to the one -degree -of -freedom response for this 
purpose is given in reference l4. Finally, it was found simple and 
expedient to obtain good estimates of the characteristic frequency Nq 

directly from the flight records of acceleration for the various airplanes 
in accordance with the method given in reference l4. 
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Variation With Type of Operation 

In reference the foregoing techniques were applied to derive 
estimates of the distribution of root-mean- square gust velocity 

from data for eight transport operations. As might be expected_, little 
difference was noted in this distribution for similar types of operation. 

A distinctive difference was noted^ however, between the results for a 
feeder-line low-altitude operation and the normal higher altitude trans- 
port operations, as might be expected. The results obtained for the 
feeder-line operation and one of the more normal transport operations are 
compared in figure 6 in order to illustrate the type of variations in the 
distribution. The feeder-line operation consisted of short-haul low- 
altitude flights over the rough terrain in the western part of the United 
States and was known to spend much time in rough air. The other operation 
was a four-engine transcontinental operation mostly flown at moderate 
altitude and was known to be a relatively smooth operation. In each case, 
the results are shown for a range of values of Uy from 0 to 20 feet per 

second. Since the individual gust peaks in rough air of a given root- 
mean-square value can be three or four times the root -mean- square value, 
the results shown cover peak gust velocities up to perhaps 80 feet per 
second. 

The curves both indicate a rapid decrease in the probability with 
increasing root-mean-square value. It is clear that the transcontinental 
operation includes more time at very low values of less time at the 

moderate values of dy, and perhaps more time at the very large values 
of dy. However, since the proportion of flight time between given root- 

mean-square values is given by areas under the curves, the significant 
areas are distorted by the logarithmic plot and direct comparison of the 
curves is made difficult. For this reason, the significant points obtained 
from these curves are summarized in the table on the upper right-hand 
side of figure 6. The table gives the percent flight time for the two 
operations in smooth air ( dy < 1 foot per second) , and in moderate to 

severe rough air (dy > 5 feet per second). The table clearly indicates 

that the transcontinental operation is a much smoother operation with 
almost twice the percentage time of the feeder-line operation in smooth 
air and only about one-eighth the percentage time in moderate to severe 
turbulence. At the very high values of root-mean-square gust velocity, 
there is some indication that the transcontinental operation may encounter 
very severe turbulence, associated with the more violent thunderstorms, 
more frequently than the feeder-line operation. The time above 20 feet 
per second is small, however, being only about l/lOO, 000 of the total 
flight time. 

The description of the gust experience in this fonm is directly 
applicable to load calculations for Other airplanes in similar types of 
operation by reversing the procedures used in obtaining these results. 
However, direct application of these results would only apply to opera- 
tions having the same type of flight plan. In order to obtain results 
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'h Via. 'h 8 ir 0 niO 2!*0 fl 0 xiTDl 0 3.nd. SLpplxc 8 Lbl 0 "to 3.3rTDi"fcr3*x*y fliglrfc plQ.ns^ il 
would 10 d 0 Sirabl 0 to dotormino the variations in thoso distributions 
with altitudo, woathor condition^ and porhaps goography. Efforts in 
this diroction aro boing mads. 


Variation With Altitudo 

In ordor to arrivo at some rough ostimatos of th 0 variation of f 

with altitudoj us 0 was mads of ths summary of gust statistics givsn in 
rsfsrsncs 5. Figurs 6 of rsfsrsncs 5 prsssnts sstimatss of ths avsrags 
gust sxpsrisncs at various altitudss that ars rsprsssntativs of contsm- 
porary typss of transport opsrations. In ordsr to sstimats ths associatsd 
distributions of root-msan-sq.uars gust vslocity^ thsss rssults which ars 
in tsrms of dsrivsd gust vslocitiss wsrs, for convsnisncs, first convsrtsd 
to accslsrations by using ths charactsri sties of a representative transport 
airplane. The charts of reference ik were then used to estimate the 
appropriate distribution foimi and scale value for the root -mean- square 
gust velocity. It was found that, for each altitude bracket, the number 
of peak accelerations for the various altitude brackets could be approx- 
imated by the following distributions for f(cJu)* 


f(au) = 0.99 


1.56 


-au/1-56 

e ' + 0.01 


-au/5.18 


3.18 


for the altitude bracket of 0 to 10,000 feet. 



1 

2 ( 0 . 34)2 


for the altitude bracket 10,000 to 30,000 feet, and 


f (au) 


1 “^/ 0*51 

2 ( 0 . 31 )^ 


for the altitude bracket 30,000 to 50,000 feet. These ^stributlons are 
shown in figure 7 . These distributions are a minor refinement of results 
presented earlier. The two -term approximation for the lowest altitude 
bracket was found desirable in order to represent the data over the whole 
range of gust velocity adequately. 

The results in figure 7 indicate that, except for the lowest 
altitude bracket, little difference exists in the distribution of Cy 
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for the various altitude brackets. Perhaps the most important points, 
which are obscured by the logarithmic scale, are the relatively large 
amount of time spent in essentially smooth air (cru < 1 foot per second) 
at the higher altitudes and the relatively large amount of time spent in 
moderate to severe turbulence (cfu > 5 feet per second) at the lowest 

altitude bracket. The time spent above 5 feet per second for the lowest 
altitude bracket is roughly five times as great as that for the higher 
altitude brackets . 

It should be remembered that these results are in terms of true 
gust velocity. If eq_ui valent gust velocities which are more directly 
related to the airplane response are used, the decrease in turbulence 
intensity with altitude would, of course, be even more pronounced than 
is indicated here . 

As a check on the consistency of these results, the acceleration 
histories were calculated for two hypothetical operations: a moderate- 

altitude transport operation and a high-altitude operation perhaps rep- 
resentative of Jet transport operations. In the calculation of the value 
of A, theoretical results for a single degree of freedom (vertical motion 
only) (ref. 9) vere used. The results obtained in this manner were 
roughly equivalent to those obtained by discrete gust calculations which 
would be expected for this single -degree -of -freedom case. 

The method of application of these results to the calculation of 
load histories for arbitrary flight plans is parallel to that used for 
the discrete gust case described in reference 5* The sequence of steps 
is: (l) The appropriate distribution of f(oxj) is selected for each 

altitude bracket. (2) The values of A are determined for each signi- 
ficant portion of the flight plan in accordance with equation (8). 

(3) Each of the distributions of T(o'u) is transformed in order to obtain 

the associated distributions of acceleration i’(cfa) relation 

f(%) = = f(ou) 

A 


where 

(4) The values of Nq are determined in accordance with equation (4a) 
or directly from flight records as described in reference l4. It should 
be noted that the value of Nq is infinite for the one-degree-of-freedom 
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case treated in reference 9 therefore the results obtained in that 
analysis cannot be used for this purpose without some modification. 

For a rigid airplane, the value of Nq appears to be larger but of the 
same order as the period of the longitudinal short -period oscillation. 

(5) The distributions of ba and the values of Nq are then used in 
equation (6) to derive the mmber of peak accelerations per second or 
per mile for each condition. (These calculations can be facilitated 
by the use of charts such as given in ref. l4.) (6) The results obtained 

in step 9 are then weighted in accordance with the flight distance in 
each condition and then summed for all conditions in order to obtain 
the overall acceleration history. 


Variation With Weather Condition 

In the preceding discussion, the variation in the gust experience 
with flight altitude was considered. Another breakdown of the gust .expe- 
rience which may be useful in some problems is the variation in gust expe- 
rience with weather conditions. Figure 8 shows estimates of the variations 
in f that have been obtained for several types of turbulent weather 

conditions. The curve labeled ”clear-air turbulence" was based on data 
obtained in flight through clear and turbulent air at the lower altitudes. 
The curve labeled "cumulus clouds" was based on data obtained in flight 
under moderate convective conditions such as represented by bulging cumu- 
lus clouds. Finally, the curve labeled "thunderstorms" was based on data 
obtained in flight in the immediate vicinity of and within severe thimder- 
storms . The distributions obtained appear to have the same general form 
but are seen to vary widely in intensity. A simple measure of the relative 
intensity of the turbulence for these conditions may be obtained by com- 
paring the root-mean-square values at the lower probability levels. It 
will be noted that, at these probability levels, the values for the cumulus 
and thunderstorm conditions are roughly two and three times the values 
of Gy for the clear-air condition. 

It has been estimated that contemporary transport operations spend 
about 10 percent of their flight time in this clear-air turbulence con- 
dition, 1 percent in the cimaulus condition, and perhaps 0.05 percent 
in thunderstorms. These results may therefore find some application in 
evaluating the effects on the overall load experience of operational pro- 
cedures which would tend to modify this weather experience. For example, 
the Introduction of airborne radar for weather avoidance may be expected 
to reduce the exposure time to the severe turbulence conditions of thunder- 
storms. Also, high rates of climb and descent through the lower and more 
turbulent altitude layers may cause a drastic reduction in the 10 percent 
of the flight time attributed to clear-air turbulence conditions and 
thereby cause a marked reduction in the overall gust experience. 
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CONCLUDING REMARKS 


The foregoing analysis has served to demonstrate that the gust 
statistics may, under reasonable assumptions, be converted into a form 
appropriate for spectral-type calculations. The significant and funda- 
mental quantity, for this pui^ose, appears to be the probability dis- 
tribution of the root -me an- square gust velocity. The results obtained 
in defining the variations of this function with type of operation, 
flight altitude, and weather condition provide at least a starting basis 
for their application to response calculations in arbitrary operations. 
These results should serve to supplement the discrete-gust techniques 
in current use and be particularly appropriate in problems requiring a 
more detailed accounting for the airplane dynamics than is possible by 
discrete-gust techniques. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics , 
Langley Field, Va., May 10, 1955* 
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Figure 1.- Sujnmary of airplane measurements of the power spectrum of 

atmospheric turbulence. 
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Figure 2.- Measurements of the power spectrum of atmospheric ti^bulence 
obtained at 300 feet from a meteorological tower (ref. 16). 
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Figure 3-- Analytic representation of the spectrum of atmospheric turbulence. 
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Figure 4.- Relations between counts of acceleration peaks and root-mean- 
square acceleration values for the homogeneous and composite cases. 
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Figure 5.- Relations between counts of acceleration peaks and the 
probability distribution of the root-mean- square gust velocity. 



Figure 6 .- The probability distribution of the root-mean- square gust 
velocity for a transcontinental and a local feeder-line operation. 
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Figure 



- The variation of the probability distribution of the root' 
me an- square gust velocity with altitude. 



The variation of the probability distribution of the root- 
j]]g aji— s quare gust velocity with weather condition. 
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